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Abstract: 

* 

This report presents the results of Part I of the dynamics 
and control analysis for the MPTS antenna and collector Interactions, 

The objectives of this part of the study have been to establish the 
basic dynamic properties and performance characteristics of the antenna 
so that the results can be used for developing criteria, requirements , 
and constraints for the control and structure design. Specifically, 
the vibrational properties, the surface deformation, and the corresponding 
scan loss under the Influence of disturbances have been studied. 
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SECTION I 


INTRODUCTION AND SUMMARY 

This report presents the results of Part I of a two-part study 
of the dynamics and control analysis of the MPTS* Antenna and Its Inter- 
action with the solar collector. The effort reported here deals with the 
dynamics problem of the antenna/collector, especially the surface deforma- 
tion and the effects of the power loss due to the warping of the antenna 
surface. The second part of this work deals with the attitude and pointing 
control problems which will be summarized In Volume II of this report. 

Although the main emphasis of this report is on the MPTS antenna, the over- 
all results of this study apply to the SPS configuration formed by a 
collector structure plus two end-mounted antennas (see Pig, I, in Section 2.2), 

1.1 GENERAL BACKGROUND 

The SPS Is the largest space system conceived to date that 
appears feasible with reasonable extensions of existing control technology. 

It represents a class of large platform-like structures that are several, 
orders of magnitude larger than any of the other large space systems 
(multiple-payload platforms, parabolic reflectors, etc.) currently in 
planning within NASA. The SPS has in common with all large space systems 
many control problems that are widely recognized within the controls 
community. These problems include attitude errors due to disturbances, 
potential instabilities due to truncated modes and other model errors, lack 
of damping. Inaccurate preflight knowledge of the vehicle dynamics, and the 
parameter variations while the system is in operation. The qualitative 
nature of these problems (model errors, concentrated stresses due to large 
actuator size, etc.) has emerged as a result of studies In the general 
area of control of large space structures. However, there is a need at 
this time to Investigate the dynamics and control problems specifically 
related to the Satellite Power System to assess performance of selected 
control concepts, and to identify and initiate development of advanced 
control technology that could enhance feasibility and performance of the 


* 

MPTS stands for Microwave Power Transmission System. 
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SPS systitm. two of ttm areas Chat have been under Investigation are 
the dynamics and control of the solar collector with the MPTS antennae 
treated as point masses and the dynamics and Ofaitrol of the MPTS antenna 
with the solar collector treated as a dynamic disturbance source* This 
report covers the most basic problems of the latter, that Is, the dynamic 
properties and Mie performance characteristics of the MPTS antenna. 

It Is a knovm fact that the losses accrued at the later 
stages of a process are more costly than those accrued at the early stages. 
This is also true In the context of the power collectlon-converslon- 
transmlssion process of the SPS system. Since the antenna power trans- 
mission constitutes the last part of the efficiency chain for the In-orblt 
operation, a great deal of emphasis has been made directly to the pointing 
accuracy. However, due to the high weight penalty at synchronous altitude 
and the huge array size (1000-meter diameter) , the structure cannot be 
made arbitrarily rigid, and the structure stiffness and hence the surface 
deformation of the anteren^ plays an Important role In the determination 
of the performance of the antenna. It Is this latter subject area which 
this study Is focussed on, 

1.2 PERFORMANCE REQUIREMENT GUIDELINE 

In References 2 and 3, the requirement of mechanical pointing 
and alignment accuracy of 2 to 3 arc minutes has *^)een considered. In 
this report, the performance requirement of 3 arc minutes for the surface 
flatness Is used as a base for discussion. 

1.3 OBJECTIVES AND SCOPE 

The combination of flexibility, huge dimensions, high distur- 
bances, and the nature of the operation make the MPTS/SPS antenna uniquely 
different from any antenna or spacecraft that the aerospace community has 
ever designed. The objectives of this part of the study have been to 
establish the basic dynamic properties and performance characteristics 
of the antenna so that this knowledge can be utilized to form a new base, 
requirement, and constraints to aid the control system design. Specifically, 
the following subject areas have been addressed based on the most up-to-date 
MPTS data. 
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It tha vtbKational pcopatclos «uch as Che mode shapes 
and natucal frequencies of Che anCenna and Cheir 
effecc Co Che SPO sCrucCure propeccies as a whole# 

The surface deformaclon of Che ancenno sCrucCure and ICs 
effecC on Che scan losses, 

3. fhe effects of disturbances Including thermal distorclons 
of the solar collector to Che flatness of the antenna 
surface. 

, Application of the linear analysis technique to extend 
Che results obtained through simulation. 

The approach employed here has been one of the time domain 
analysis techniques, i.e,, a combined modeling, analysis, and computer 
simulation approach. 

In Section 2, the necessary mathematical formulations of the 
model and the performance related terms and quantities are presented. The 
gti'uctural flexibility properties, the estimates of the disturbances, the 
thermal effects, the dynamic Interactions, and the main results of the 
simulation and their linear extensions are presented in Section 3. 

1.4 STUDY RESULTS AND CONCLUSIONS 

The main contribution of this part of the study is that It has 
established the performance characteristics of the antenna, especially the 
local warping and its effects on the antenna scan losses. The performance 
characteristics are plotted as time histories of a number of quantities and 
they have also been summarized in tabulated forms in Section 3. In 
addition, these results have been extended In the parameter space as shown 
in Figs. 23, 24, and 25 (these figures are reproduced here for quick 
reference) . 

Throughout this study, a great deal of insight and experience ^ 

have been gained such as how certain types of vibration modes react to 

* 

given types of signals (forcing functions) and their position of execution; 
how modal dominance varies with the frequencies and the properties of 
the forcing functions; the importance of signal shaping and timing; the 
distribution of modal energy and its time dependence, etc. All 
of these are invaluable information for the control system design, 
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Computer simulation programs have been developed which have 
proved to be an effective and necessary tool for this study and it will 
be an essential part of the facility for our future work. These programs 
are listed in Appendices A and B for reference purposes. 

The major conclusions of this study are as follows: 

1, The structural damping is essential for the maintenance 

of a flat radiation surface. This is especially Important 
for dissipating the otherwise prolonged presence of high 
surface transient caused by disturbances. The decay rate 
of the surface deformation in general agrees with that of 
the amplitude of the dominant mode. A structure damping 
ratio of .005 has been used; a higher value will be more 
desirable, 

2. The most significant sources of disturbance that affect 
the antenna performance perhaps come from the dynamic and 
control interactions between the antenna and the collector. 

The effects of these disturbances may be reduced Blgnlflcantly 
by 

a. decoupling the antenna from the collector motions, 
especially the translational motions; 

b, actively controlling the collector bending motions 
at the interface boundary. Large amplitude 
oscillations at the collector tips must be actively 
suppressed. 

Figures 23 and 24 show the performance characteristics of 
the antenna surface as a function of the amplitude of 
oscillation (cosine function) for collector bending and 
torsional motions, respectively. A suddenly applied 
bending oscillation of slightly over 1/4000 of the 
length of a 20 km collector will cause a maximum RMS 
local slope angle of 0.05® or 3 arc minutes. 
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External disturbances such as the gravity gradient, solar 
pressure, etc., will not have a direct Impact on the surface 
flatness with the possible exception of the thermal dlator- 
tlon of the collector structure. High temperature gradient 
with short thermal lag time after shadowing, can set the 
collector to sudden bending oscillation, Unless low 
coefficient of thermal expansion (GTE) graphite composite 
material Is used, this could be a source of serious problems. 

^Figure 25 shows the estimates of the collector bending 
amplitude as a function of the GTE for a number of tempera- 
ture gradient values. By using this estimate with the 
performance data of Fig. 23, a suitable GTE value may be 
determined . 

Signal shaping and timing are of critical importance to the 
performance of the MPTS/SPS system. The shape of an 
actuator force should be designed so that the least amount 
of energy may be absorbed by the flexible modes. The cut 
off time of a thrusting force determines the amount of 
energy left in the flexible modes. To reduce the level of 
surface vibration, the modal state of the dominant mode 
should be monitored so that the cutoff of a major thruster 
should be timed such that the modal energy of this mode Is 
at Its minimum. 

The MPTS/SPS has the following structural properties; 

a. The natural frequencies of the antenna and the 
SPS are 2 to 3 orders of magnitude greater than 
the orbital frequency. That means that these 
modes and orbit will not couple, 

b. The frequency of the lowest flexible mode of the 
antenna Is about one order of magnitude greater 
than that of the collector. However, the first 
3 flexible modal frequencies of the MPTS overlap 
with the 8th of the collector, arid the excitation 
of the latter should be avoided. 
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c» Due to Its location, the antenna mass will have more 
Influence on the SPS normal frequencies than the same 
mass distributed over the solar collector. The antenna 
mass Is about 30% of the total SPS mass, its existence 
has decreased the fundamental frequency by a factor of 
0.4 and Its effect on the higher frequencies la less 
[5] , The stiffness of the coupling structure also 
affects the normal frequencies. The significance of 
these effects to design alterations are yet to be 
decermined . 

d. The most significant factors that affect tlie modal 
properties are the geometrical parameters of the 
structure. 


AMPLITUDE OF ACCELERATION FORCE, N 
1367 2733 4101 5468 6834 8201 





S. 



Figure 23, MPTS Surface Deformation vs. Amplitude of 
Collector Bending Oscillation 
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Figure 25 . Co! 






SECTION 2 


DYNAMIC MODEL AND PERFORMANCE RELATED EQUATIONS 


2.1 TIIS DYNAMIC EQUATIONS AND ENERGY OF CONCENTRATED MASSES 

Dynamic Equations 

In Reference 1, finite eJement: niodels of Che MPTS antenna 
have been discussed. In this sectltnii only the key elements of the model 
that have been Implemented In the simulation programs are discussed. 

The model consists of 20 decoupled second order differential 
equations representing the modal amplitudes of the 6 rigid body modes 
and the 14 lower order flexible modes, 

Mq + Kq “ ^"^f (1) 

where q Is the 20-^modal amplitude vector; f, the 498-force vector; M and 
K are the 20x20 generalized mass and generalized stiffness matrices (both 
diagonal), respectively; and Is the matrix of eigenvectors of dimension 
498x20, For convenience, (1) Is written in the following form, 

q + Aq £^f (2) 


where 






\J 


(3) 


where Is the ang'.ilar frequency of the kth mode for k*“7,. . ,,20 and 
" 0 for k«l,,..,6; and 

(M Is diagonal) (4) 


In this study, structural damping is also considered. Since there Is 
no available information, a damping term, with “ .005, has been added 
to the kth flexible mode as follows, In the expanded form, 

**k “ ^ ’ k”l,...,6 


(5a) 


(5b) 


\ 2C4ij^q + £, k»7 , , . . >20 


where <jij^ la the kth column of ♦ and is the kth diagonal element of 

Since these equations are linear time invariant, analytical solutions 
are readily known. Compact solutions were therefore used in the simulation 
programs Instead of the more time consuming methods of numerical integrations. 
In the following, let T be the interval of each computation step, and 
♦ - ■■ cosiji, “ sinij), and q(n) “ q(nT), then in each 
computation step, the following is evaluated, 


For k “ 1, . . . ,6 


q(n+i) 


1 T 


q(n) 


T^/2^ 


m 




+ 


q(»'+i)j 


0 1 




T ^ 




For k “ 7 , . . . , 20 
q(n+l)" 


q(n+l) 


® e 


-C^nT 




w. 




q(n) 

q(n) 


I - (c^ + V' 

\ ^ 


£(n) (6b) 


where f(n) is the force vector evaluated at t “ nT and n ■ 0,1, . , . ,tj/Tj 


t^ is the final time. 


With the availability of qj^, the displacement of node i, 
i “ 1,... ,166, due to elastic deformation alone, at t = nT, can be 
evaluated as the linear combination of q^^, i.e.. 


(7a) 


20 

“l<”> - ♦kl 

20 

''!<"> ■ ♦k(l66tt) 1k<"> 

20 

•t<“> ■ ♦k(332+l) ’k<"> 

and the displacement due to the rigid body inodes, 

6 

"i<“) ■ ♦ki \<"> 

^i<“> ■ Ji ♦kaee+i) ’k'") 

■ X *>‘(332+t) ’k<"> 

V 

The actual position of node i, at t ■* nT, is 
Xi<n) - + u^(n) + Uj^(n) 

y^(n) “ v^(n) + u^(n) 

z^(n) « + w^(n) + Wj_(n) 

where Y^, are the coordinates of node 1 when the system is at rest. 
Equation (9) was used to generate the antenna surface and the local slopes 
with the rigid body motion suppressed, i.e.. 


Xj^(n) 

“ X^ + u^(n) 

(10a) 

Yi(n) 

« Y^ + v^(n) 

(10b) 

\{n) 

“ Zf + w^(n) 

(10c) 


where (Xj^(n), y^^Cn) , z^(n)) may be interpreted as the coordinates of node i 
at t “ nT projected onto the body frame, whe?."eas (x^(n) , y^^(n) , z^(n)) is 
that expressed in terms of an inertia frame. 

Energy of Concentrated Masses 

The body energy of concentrated masses is of great interest to 
this analysis Since it is a good measure of modal activities of the antenna 


(7b) 

(7c) 

(8a) 

(8b) 

(8c) 

(9a) 

(9b) 

(9c) 
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undei: disturbance. Let E be the sum of kinetic and potential energy of 
the antenna • Recall that the antenna Is modelled as a collection of 
concentrated masses which are connected and having finite stiffness. E 
may be expressed 

E » j X 'mx -f ^ x^kx (11) 

T T 

Since X " and !s£* becomes, 

E - I Mq + i q^' Kq (12) 

Recall that M and K are diagonal. E may be written In terms of its modal 
components, 

K-1,,,.,6 (13a) 


\ ■ I + “k^ ’k^) - 


7,..., 20 


(13b) 


Let Ej^ and E^ be the rigid body energy and the flexible body energy, 
respectively , then 


and 


E « i y M, 4 ^ 

1 „ /• 2 ^ 2 2 . 

k«7 


E =. + Ej 


(14a) 

(14b) 

(14c) 


Since qj^ and qj^ are computed In each step, the energy Ej^, 

E , E_, and E can be readily evaluated. It will be Illustrated later 
R i* 

that the modal energy is Indeed a good measure for determining the relative 
dominance of the modal activity. Further, it may be used as a measure for 
signal shaping and thrust on-off timing. 


2.2 LOCAL SLOPES OF THE SECONDARY STRUCTURE 

The basic architecture of the MPTS antenna consists of two 
parts, the primary structure and the secondary structure (Fig. 1). The 
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primary scrucCute Is a tatrahadral planar truss structure which serves 
as the base support to the system and provides the required structural 
stiffness. The secondary structure Is a "finer” structure which Is mounted 
on top of the primary and It provides a base for the Installation of 
subarrays of the antenna. Due to the huge size (1000 m dla« x 130 m depth 
- nominal) and the weight penalty of a satellite to be built and trans- 
ported to a synchronous orbit, the stiffness of such on antenna cannot be 
made arbitrarily high. A series of questions may be asked, such as, 

(1) how stiff the structure will have to be In order to meet the operational 
requirement?; (2) Is the configuration under study stiff enough?; (3) Is 
there any trade off between structural stiffness and active shape control? 
Before these questions can be addressed, the relationship between the 
structural deformation and the ff.>tenna transmission loss must be established, 
This Is the subject of the next subsection. In this subsection, the local 
slopes of the antenna radiation surface Is derived. 

For structures of finite stiffness, surface deformation will 
occur if disturbances or forces are applied to the structure. Accompanied 
with the local displacement is the change of the direction of the normal 
to the local surface, or the local slope. Since the antenna Is retro- 
directive phase conjugate, the transmission efficiency lo not sensitive 
to the variation of path length or the surface dispositions but it Is 
very sensitive to the local slope changes. Therefore, It Is Important 
to evaluate the local slope variations. 

Since the finite element model being used does not provide the 
details at the subarray level, one can only compute the slopes of the 
secondary structure. In fact, the secondary structure is not modelled 
either; its slopes are computed by assuming that the local surfaces are 
flat and each one a part of one of the 61 plane surfaces. Referring to 
Fig. 1, the secondary structure consists of 61 hexagonal substructures 
referred to here as the hexagonal planes. Each hexagon plane is supported 
at the three symmetrical vertices by the supporting pins on the primary 
structure such that there is no direct contact between the primary and 
the secondary structures and that there is no contact between the secondary 
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Substtuctuces themeelveo* This Is done to teduce the Interactions 
between structural components. It may also be Important to point 
out that the grid points assigned In the model coincide with the 
supporting pins (actually each grid point represents the three 
closest pins and Is positioned at the geometrical center of the three 
pins (refer to Figs. 2 and 3)). 

With this background in mind, the local slope angles can be 
readily derived. There are at least two ways for computing the slope 
angles, ( 1 ) using the equation of a plane, and ( 2 ) using vector product. 

Let (xj^, yjL* z^). (X 2 » 5^2 * ^2^» ^^3* *3^ 

coordinates of the three grid points associated with a typical hexagon plane 
1 (the index 1 and the grid point index j were dropped) . In the following 
the bat in ”7'' that signifies the grid point position with the rigid body 
motion suppressed Is also dropped. 

Method 1 - Plane Equation Approach: 

Let a, b, and c be the constant coefficients, and the equation of 
a typical plane may be written as, 

ax + by + cz » 1 (15) 

The coefficients may be determined by substituting the three points Into 
(15) and solving the three linear simultaneous equations, or 


a 


'^1 ^1 "l' 

-1 

* 1 ' 

b 

' as 

^2 ^2 ^2 


1 

c 


X 3 Z 3 _ 


1 


To avoid numerical Iterations , compact expressions for the matrix 
inversion may be used. Since the coordinates of the three points are com- 
puted at each step, the values of a, b, and c can be evaluated accordingly 
and the slope angle error, 0 , is readily obtained. 


(a) 75 Grid Polnca on Array Side of Primary Structure (Only 13 with 

ID shown) . 


SUFPORTING PINS 



(b) Secondary Structure Supporting Pins. 


Figure 2, Grid Points of the Finite Element Model and the 
Secondary Structure Supporting Pins 
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Note Shat the matrix inverse in (16) does exist , since, in our application 
the three points ore olways non*>collnear, However, this roett od could lead 
to computational difficulties. Referring to Fig. 4, 1/a, 1/b, and 1/c 
arc the three intercepts of the plane on the X», Y-, and Z-axis, 
respectively, hence, as the plane becoming parallel to the XY -plane, 
a ^ 0, b *»• 0 end c >*■ «, Since 0 is a very small angle (and time varying) , 
double precision computation will be required for the evaluation of (17) , 

Method 2 - Vector Product Approach! 


This approoch is simpler than the first method and it requires 
fewer computation steps. Let and A^ be defined as follows. 


A 


I 


X2 - Xj_ 
^2 - ^1 
. ^2 “ h 


(18a) 



(I8b) 


Let their vector product be (A B C)^, then 


where 


A 

B 

C 




A. A» 
1 2 



the slope angle is 


0 => cos 


-1 






B*’ + C 


(19) 


( 20 ) 
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For the same reason discussed in Method 1, double precision computation 
will also be required in evaluating (20). Method 2 was used in our 
computation. 

The slope angle, 6^(n), is a function, and its value varies 
with time (n) and space (1) . The individual values of the slope angle 
will not mean much for the antenna as a whole. Other performance measure- 
ments must be defined. Three measurements are used here. 

The Maximum Slope Error 

e^(n) - Max 0^(n) (21) 

ie{l,. . . ,61} 

6^(n) is the greatest slope angle of the 61 hexagonal sub- 
structures, at time t = nT; therefore, it represents the worst case. 

The RMS Spatial Average 



It I (e^(n))2 

1=*1 

0^s('^) fluctuates with time which is a good 
at a typical time t *= nT. 

The RMS Time Average 

/sir ra=l 

0 ,j„(n) is a smoother measurement especially when n becomes large. This 
Krlb 

measurement is useful when one asks how is the antenna doing so far. 0„,,„(n) 

RMS 

becomes less and less sensitive to the current sample as n increases, 
which is a common defect of the averages of this kind. This problem may 
be solved by applying a window to the samples, i.e. , one only uses the N 
most recent samples in the average. Another way to increase the sensiti- 
vity is to assign a heavier weight to the most current sample, for instance, 
assign 10% weight to the most current sample and 90% to the rest as a 
whole when n ^ 10. This latter approach does not require additional 
storage space. 


( 22 ) 

measurement of the antenna 
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2,3 THE MICROWAVE TRANSMISSION SYSTEM SCAN LOSSES 

Since the design of the MPTS antenna Is still In the evolving 
stage, the data and assumptions made In tills section are far from final. 
Nevertheless, the results of this and the subsequent sections will still 
be valuable Information for the purpose of developing performance require- 
ments. 


The MPTS antenna Is a phase array antenna which consists of 

a large number (say, on the order of 10,000) of array elements. Array 

elements are Installed on the 61 secondary substructures. Each element 

2 

has an aperture area of about 108 m , or for a rectangular aperture the 
aperture size Is about 10,39 m. 


The array Is active rctrodlrectlve so that it Is made Insensi- 
tive to the path length variations but It is quite sensitive to angular 
deviations from the llne-of-slght . The latter requires accurate pointing. 
Pointing will be achieved by two control systems, the mechanical pointing 
control which Is required to stay within a few arc minutes of the target — 
the rectanna, and the fine pointing control-electronic beam steering which 
Is required to steer the beam to within several arc seconds. 


Let D be the aperture, 0 be the slope angle (or angular pointing 
error), and X be the wavelength of the microwave power. Since the frequency 
Is 2.45 GHz, the corresponding X Is .1224 m. For uniformly Illuminated 
antenna with rectangular aperture, the radiation pattern is 


sin (J^ sinO) 


ttD 


alnO 


. sin (^) 


(24) 


the approximation is valid for small 6 . Equation (24) prescribes the 

patterns for both the array and the elements; for the array, one substitutes 

D„ “ 1000 for D and for the elements, use D„ ■ 10.39 for D. From this one 
a r e 

can see Immediately that the array pattern is much narrower than that of 
the element. More precisely, the half-power-width (the beam width at the 
half power points) of the array is about 100 times narrower than that of 
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t\\a Qlcment. At the hal£»powec~p(>lnt irDO/X >* 1.3916 rad and the correa- 
ponding and 0^ are .0031** and .299", respectively. Figure 5 shows the 
antenna and element patterns. 

* 

In the following it is assumed that the electronic beam 
steering system points the beam at the rectenna perfectly and the 
discussion is concentrated on the effects due to surface deformation of 
the structure, i.e., the effects of scan angle variations on the antenna 
efficiency. 


Scan Loss 


The amplitude of the array pattern varies with the scan aiitgle. 
The field strength decreases rapidly as the scan angle increasingly 
deviates from the normal of the array surface, The governing relationship 
is that the array beam amplitude is prescribed by an envelope which is 
determined by the contributing element patterns. In the case that the scan 
angles are the same for all the subarray elements, the envelope becomes 
the element pattern itself. This is illustrated in Fig. 6. The loss of 
power due to scan is called the scan loss. 


Let be the total power, the maximum power of element i, 
and sin corresponding element pattern, where the argument 

G O X 


sin X. 2 

- I P. (r~^^ 

^ i ^ ’^i 


(25) 


a. Rigid array surface 


In this case all x^ ■ x or * 0, 
sinx^2 ^ 
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where is the maximum power of the array. The scan 
loss Pj^ and Che percent loss are, 


P » P [1 - 

^L TM • X '' J ’ 

P.% - [1 - X 100 
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Figure 5. Array Pattern and Element Pattern for Rectangular Aperture 
Array with Uniform Illumination 



Figure 6. Beam Amplitude vs Scan Angle 


The percent scan loss Is plotted In Fig* 7 which 
Illustrates how critical the scan angle is to the power 
output. 


h. Uniform Illumination 

For uniformly illuminated array, - P, 
„ slnx. „ P„ ^ slnx. „ 

■ p I - (-fi) I (-^)' 


where N is the total number of elements and 
^ U sinxj 

V X - [1 - ^ I (~- 
^ ” k-l ’‘i 


-i)^] X 100 


(27) 


c. The RMS approximation 

Assume that the root-mean-square value of the scan angle, 

6 is known, where 
RMb 


l^MS 


/ITv 

i"l 


The corresponding Xj^^g can be easily evaluated: 
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then 
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and the percent scan loss is, 
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JJote that in general F,j, end hence P 
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However, our simulation shows that for many test cases 
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Figure 7 • Scan Loss of a Rigid Array 
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they are quite close. Therefore, may bo used 

Instead of P^%. The former requites much fewer computation 
steps. Note that the 6^^ here Is the RMS spatial average. 

d. Gaussian Illumination 

In order to reduce the side lobs of the array pattern, 
the elements may be Illuminated non-unlformly . One way 
to accomplish this Is to tailor the Illumination of 
the element according to a Gaussian curve. In this case, 

p . p “ (31) 

o 

The parameter a may be determined by the taper at the edge 
of the array. For a 10% taper, l.e., P/P^ ■ 0.1 when 
R * R^, a •» 2.3, A 10-step taper that approximates a 
Gaussian curve was proposed In an earlier study [2]. 

From (25), P^ becomes, 

K “ I P. (32) 


„ slnx 

p I ( — ^ 
° 1 *1 


R. 2 


where R^ Is rhe distance from the center of the array to 
the element, and R^ the radius of the array aperture. 

Of the last three cases, only the first two were Implemented 
(refer to Appendix A) . 
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SECTION 3 


SURFACE DEFORMATION AND ITS EFFECTS ON THE ANTENNA EFFICIENCY 
DUE TO DISTURBANCES AND DYNAMIC INTERACTIONS 

This section addresses the problem o£ surface deformation 
and Its effect on the microwave power transmission efficiency In quanti- 
tative terms. Specifically, the question of how significantly the disturbances 
will impact on the antenna *s flatness and how much disturbance the antenna 
can stand without Incurring excessive power loss will be answered. The 
results were obtained through extensive computer simulation, 

Before the main results are discussed, the modal characteristics 
of the antenna and how these affect the overall system modes Is briefly 
described and then the magnitude of disturbance and dynamic interaction 
forces are estimated which form the basis for the magnitude of simulation 
Input . 


3.1 THE MODAL CHARACTERISTICS OF THE ANTENNA AND THEIR EFFECTS ON 

THE OVERALL STRUCTURE OF SPS 

The fundamental properties of a large space structure such as 
the MPTS antenna that have significant effects on the control system design 
are the structural vibrational properties which are governed by the funda- 
mental frequencies and the mode shapes. The former affects the controller 
bandwidth design and the latter affects the sensor/actuator location 
selections . 

The first ?0 natural frequencies of the structure are summarized 
In Table 1. Of the 20 modes, the first 6 are the rigid body or zero 
frequency modes and the rest are flexible modes. Figure 8 shows the 
mode shapes of the first 14 flexible modes. The first two modes in Fig. 8 
exhibit astigmatism bending, the third one Is a defocus mode, and the 
fourth mode (mode 10) exhibits trefoil bending characteristics. The rest 
of the flexible modes have either the similar modal characteristics of 
the first four but with higher frequency patterns or a combination of 
them. 
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The mode ihepes of the first: four flexible modes are in 

general agreement with that discussed in Ref, 3 where only the first 

four flexible modes were illustrated. However, the modal frequencies 

are about 4 to 7 times greater in Ref, 3 than those listed in Table 1. 

It is also noted that the nominal mass used for the antenna in our model 
6 6 

is 15 X 10 kg whereas it was 8,58 x 10 kg in Ref, 3, a factor of 1,75. 

Since the generalized stiffness is proportional to the mass times 

2 

frequency squared (mm ), the stiffness is about 11 times greater for 
modes 7, 8, and 10, and greater for mode 9 in Ref, 3 than these in our 
model. 


Figure 9 shows a comparison of the natural frequencies of 
the MPTS antenna, solar collector, and the first five frequencies of the 
flexible modes of the SPS. The following conclusions may be drawn. 

1. The fundamental frequencies of the SPS and the antenna 
are 2 to 3 orders of magnitude greater than the orbital 
frequency. This means that the couplings between the SPS 
flexible modes and the orbit will not be significant. 

2. The fundamental frequency of the antenna is about one 
order of magnitude greater than that of the collector. 
However, the frequencies of the first three flexible 
modes overlap with that of the 8th mode of the collector 
which Indicates possible couplings between these modes. 

Due to the fact that higher frequency modes are less 
likely to be excited, the likelihood of their exciting 
the lower modes of the antenna are also reduced. 
Nevertheless, caution must be exercised in the controller/ 
actuator design so that mode 8 of the collector will not 
acquire excessive energy at any time. 

3. The antenna does not decrease the SPS frequencies signifi- 
cantly nor does the coupling stiffness [5] . This is 
because of the superior size and mass of the collector 
over that of the antenna. The implications of this 
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obsctvaclon is that tha doslgn altetatlons o£ the 
antenna and the coupling will not significantly Impact 
on the vibrational properties of the overall system. 

4. The moat significant factors that affect the modal proper- 
ties arc the geometrical properties of the structure for 
both the antenna and the collector, l.c., the size and 
the relative size (e.g,, aspect ratio) of the dimensions 
of the structure. 

The 3-D plots of the mode shapes was done by the program listed 
in Appendljc B and the required dynamic responses were computed by the 
program listed In Appendix A. 

3.2 ESTIMATES OF DISTURBANCES AND DYNAMIC INTERACTIONS 

Disturbances may be catogqrlzed as external and onboard. 
External disturbances include gravity gradient torques, solar pressure, 
magnetic torques (due to the interaction between the Earth's magnetic 
field at the synchronous orbit and the onboard current loops), etc. 

The onboard disturbances include, for instance, the microwave recoil, 
dynamic unbalance torque, current interaction torque, thruster firing, 
etc, Qualified as both external and onboard is the thermal bending distur- 
bance due to temperature gradient of the structure. Excluded from either 
category is the dynamic interaction between the array and the collector 
such as the forced motion of the antenna due to collector bending. 

The effects of some of the disturbances can be reduced by 
careful design. For Instance, the current interaction torque may be 
reduced by careful arrangement of the conductors so that the conductor 
pairs are ccplanar or by using coaxial cables, Thruster impingement may 
be reduced by using glmballed pairs and by careful selection of their 
locations, etc. Solar pressure and microwave recoil forces (and torques) 
are rather small to cause any flatness problems of the antenna. They 
may create station keeping problems more than anything else. Therefore, 
it is sufficient only to look into the effects of gravity gradient torques 
and the forces due to collector bending. 


I 


Gravity Gradient Torques 

The gravity gradient torque about the x~axls (refer to Fig* 10) 
Is, 

Tx - - f (I2 - ly) sin 2 (| - 

- - I (I2 - ly) sin 2 0^ (33) 

where x-, y», and z-axes are the body axes of the antenna, x and y are 
the Inplane axes and z-axls Is normal to the array surface. 1 and I„ are 

y ^ 

moments of inertia and their values are, 
ly - » 7.24 X 10^^ Kg-m^ 

I, « 1.45 X 10^^ Kg-ra^ 

and the orbital rate Is 7.272 x 10 rad/sec. 

Assuming that the orbit Is IncUncd at 7.5°, the satellite 
will oscillate of + 7.5° about the equatorial plane once a day (refer to 
Fig. 10). For a rectenna located near the equator, the attitude of the 
array will have to vary + 1,34° about the local vertical once a day. 

For a rectenna located at 45°N the angular deviation will vary between 
6,04° and 7.53° from the local vertical. For an equatorial orbit, the 
angular deviation will be a fixed bias of 6.83° for the 45°N rectenna. 

Of all these cases, the maximum deviation angle from tl»e local vertical 
Is 7.53°, the maximum bias angle is 6.83°, and the maximum cyclic angular 
amplitude is 1.34°, and the corresponding gravity gradient torques are, 
respectively, from (33), -1496 N-m, -1360 N-m, and -269 N-m. These torques 

are Instabilitlzlng disturbances, i.e. , 0^ ■ 0 Is not a stable equilibrium. 
In all of these cases, control effort will be required to offset 
gravity gradient effects In addition to tracking (attitude 
correction to compensate latitudinal motion due to nonzero Inclination 
and to compensate longitudinal motion due to nonzero eccentricity) . 

The control problems will be discussed In Volume II. 

The problem of array flatness Is the main concern here. It will be clear 
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ANTENNA 






Figure 10. Gravity Gradient at Three Orbit Positions oi; 
the MPTS/SPS Oscillating with a 24-hour 
Period. 



later that the warping of the array surface due to gravity gradient Is 
relatively small and Insignificant due to the relatively low disturbance 
level and the low frequency nature of this torque. 

Collector Bending Forces 

Collector bending motion may be caused by many reasons, for 
Instance, thermal effect, thruster firing for station keeping, maneuvering, 
attitude acquisition and control, etc. Unless the translational and 
angular motions of the collector at the Interface boundaries are controlled 
or minimized, they may be the most significant causes of the antenna 
warping and power loss. 


a. Collector bending due to thermal effect 

In order to estimate the magnitude of the thermal bending 
distortion, a collector configuration must be assumed. 

To be consistent with the work of Johnson Space Center [5] 
and for the purpose of comparison, the following analysis 
employs a full configuration Identified as 20D4 In Ref. 5. 

Let A and d be the length and depth of the collector 
structure, respectively, and their values are, A - 20,000 m, 
d = 400 m. Assume that the temperature Is uniform on each 
section parallel to the collector surface and the tempera- 
ture gradient Is uniform along the depth. Let AT be the 
temperatuvvf! difference between the front (solar blanket side) 
and the back surfaces of the structure. Let a be the 
coefficient of thermal expansion and A\ be the displacement 
of the antenna Interface from the boundary of the collector 
structure. Under these simplified conditions, the bending 
displacement, H, Is, referring to Fig. 11, 


H = 


d 

oAT 


(1 


[1 - COs(' 


fa AT 
2d 


•)] 


+ f * sln(- 


faAT 

2d 


•) 


(34) 


The displacement with respect to the center of mass, H^,, Is 
= .444 H (35) 
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For a temperature dl£ Terence AT * 100 “C, the correspomllng 

values of H oiul for three values of a are ahowu in 
c 

Table 2. 

I* 

The coefficient of 10“'* m/m/“C is in the lower range of 

metal. For instan\';a, pure aluminum !ms a tliermal exi)ansion 

coefficient of 2.36 x 10"^/®C at the temperature range of 

293 K to 393 K and /illghtly less for aluralmim wrougl^t alloys. 
“6 

The values of 10 and 10 are in the range of graphite 
composite material. 

The thermal disturbance on the solar collector is the worst 
at the vernal equinox and the autumnal equinox whore the SPS 
will stay in the full shadow the longest period of time. 

Tlio thermal transient occurs wlien the collector moves into 
the sluidow and again when it moves into the sun light from 
the shadow. Since the transit time during the penumbra is 
about 1.63 seconds, which is about 1/570 of the period of 
the first bending mode of the SPS structure, the effect duo 
to tlie penumbra is ignored. For a similasc argument, the 
thermal lag time is also Ignored since the lag time is 
rather small duo to the sparsely distributed structural mass. 
As a result of these idealized assiunptions, one may construct 
this scenario. The tips of the structure arc bent away from 
the sun while the structure is heated by the full sun, and 
the equilibrium condition is reached. Suddenly the collector 
enters the shadow and the temperature gradient is reduced to 
an insignificant value, and ttie distorted structure is 
released from its retonslon force and starts to bounce 
back. This causes an oscillatory motion at the tips 
about the new equilibrium state. This motion applies a 
sinusoidal cosine translational acceleration or forces at 
the antenna Interface. 














The above scenario Is drawn to help in the determination 
of the proper forcing function. The real Interest here 
Is not In the period when the SPS Is In the occultatlon, 
but Is in the period after the uccultatlon. With the 
above discussion, one can construct a similar but reversed 
scenario; that Is, the structure starts a bending motion 
after It enters full sunlight. The only difference Is 
that this time the equilibrium surface Is a bent surface 
whereas in the shadow It Is a plane. 

Let Z be the displacement of the antenna Interface, then 

Z = -Z cos 0 ) t (36) 

ms 

and 

•• 2 

Z = Wg Z^ cos Wgt (37) 

= A cos w t 
m s 

and the equivalent force applied at the antenna Interface 
will be, 

2 

F = m 0 ) Z cos 0 ) t (38) 

asm s 

where Z = H , (o = 2ir, and f = .00675 rad/sec and is the 
m s ' s g 

first bending frequency of the SPS. m = 15 x 10° kg is 

a 

the nominal mass of the antenna. Table 3 shows the values 

of Z , A , and F for the three values of a. 
m m . m 

The effects of the thermal disturbance on the antenna's 
surface deformation are examined In the next subsection. 

b. Collector bending due to other disturbances 

The collector structure bending caused by disturbances 
other than thermal effect may be significant. Since an 
accurate estimate requires the details of specific designs, 
therefore, instead of being specific about the disturbances, 
a more generic approach is taken here. Since it is not 
inconceivable to anticipate a Z^ of 1/2000 of the length 
of the collector structure, forces corresponding to a 
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range about this Z were considered in the simulation. 

m 

Both sinusoidal and step functions were applied. 

Torsional Motion of the Collector at the Interface 

The first torsional mode (see Fig, 12) of the SPS is the third 
flexible mode which has a natural frequency of .00359 Hz or .0226 rad/sec [5]. t 

For a 1® torsional oscillation at the collector's tip, the corresponding 
torque applied at the antenna interface will be, 

T ■ I 0 )^ 0 * 6.4528 X 10^ N-m 

m m 

The corresponding force exerted at the antenna interface will be 
F " 6452. 8N, Sinusoidal and step forcing functions were used in the 
simulation. Note tliat a 1“ torsional oscillation at the tip of the 
collector cotrespondB to a translational oscillation of 43 » 6 in at the 
extreme boundary (point a in Fig. 12) and 8.7 m at the antenna Interface, 
which is Intolerable during the normal operation of the antenna although 
it may occur during maneuvering operation. 

3.3 SIMULATION, RESULTS, AND DISCUSSION 

3.3.1 Types of Forcing Functions 

The performance of the array structure was tested (via computer 
simulation) for a number of signal types that approximate the disturbances 
or maneuvers it may encounter during operation. 

The forces and torques applied are step functions, rectangular 
pulses, and sinusoidal functions. The first two kinds are for maneuvers 
and the sil^usoidal forces are for simulating dynamic interactions induced 
by collector boundary oscillations. In this later case, both sine and 
cosine functions were applied; the former represents a gradually applied 
cyclic acceleration whereas the latter , a cyclic shock load . 

In all the cases , the forces were applied at six positions on 
the primary structure as shown in Fig. 13. At each of the grid points 
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6025 and 2025 the force was applied and at each of 6045, 7045, 2045 
and 3045 the force F^/2 was applied. Although the 6-polnt distributed 
case Is more realistically representative of the antenna and collector 
Interface, the results show no significant difference from 2 points 
(6025 and 2025), Torques were simulated as couples. All the forces 
considered are normal to the antenna surfacci therefore, the results 
represent the worst case. Note that the antenna/collector Interaction 
will have the greatest effect on the antenna surface when the normal 
of the two surfaces are In parallel. 

3.3.2 Simulation Time Considerations 

In most of the cases the simulation time of 1000 seconds 
was used. This time covers at least one complete cycle of the lowest 
bending mode of the system (the longest modal period is 930.8 sec.). For 
the torsional acceleration, 600 seconds were simulated which covers 
more than two complete cycles of the excitation torque (the first torsional 
period is 278 sec.) 

The computation step sizes of 2,5 seconds and 5 seconds were 
used which represent about 1/20 and 1/10 of the dominant bending modal 
period, respectively, and about 1/12 and 1/6 of the dominant torsional modal 
period (refer to Table 1) , Step sizes greater than 5 seconds will be 
undesirable due to the increased modulation effect on the output. The 
improvement of 2,5 seconds over 5 seconds has not been significant. Note 
that the solution algorithm can tolerate much greater step sizes than 
the numerical Integration algorithm. For step functions, exact solutions 
will be obtained regardless of the step size used. 

3.3.3 The Simulation Outputs 

The outputs of a simulation run are plotted as time histories 
of all the Interested quantities. Tlie plot consists of the modal amplitude 
of all the flexible modes , the kinetic and potential energy of the indivi- 
dual flexible modes, the sum of the flexible modes, the rigid body, and 
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Che sum o£ all Che modes; (Che cnecg/ ploCs are parClcularly useful for 
IdenClfylng domlnanC modes); Che maximum ouC of plane dlsplocemenC; Che 
Chree kinds of local slope angles, l.e., Che rooC-mean-square spaClal 
average (RMSL), Che rooc-mean<-square running average (RMST), and Che 
maximum slope (SLOP). The antenna scan losses due Co slope devlaClons 
from Che llne>’ 0 f-slghc were computed and ploCCed In Cvo ways, the first, 
labeled as SLOSS, Che scan loss of each of the 61 surfaces compuCed 
Individually, and their effect summed; Che second, SLRMS, computed by 
using the RMS average of Che entire array. Although these are distinct 
terms, their values are practically the same. 

In addition to the Information plotted. In the output 
listings Che surface of the maximum slope and the grid point of maximum 
displacement are listed at each time step. 

3.3.4 The Surface Deformation and the Scan Loss 

A. Surface Flatness Subject to Collector Bending Osclllstlon 

The effects of collector oscillation of an amplitude of 

1/2000 of Its length, or 10 m, are discussed here. From Eq. (38), with 

6 

« 10 in, m^ “ 15 x 10 k§, « ,00675 rad/sec, the amplitude of the 

acceleration force, F will be 6832.8 N, or the forced acceleration Is 
.0465 X 10~ g. Figure 14 shows the responses of the structure. By 
examining the modal energy (Fig. 14 (e)-(h)). It Is quite clear that the 
first three flexible modes (7,8,9) are dominating. The oscillatory 
property of Fig. I4(a)-(d) shows basically the vibrations of the dominant 
modes whose frequencies range from .0189 to .0196 Hz. 

The surface deformation under this disturbance has exceeded 
the allowed limit as the RMS spatial average slope angle has exceeded 
0.09° (or 5.4 arc. min.) shortly after the force has been applied and It 
still reached about .07“ after 930 seconds as shown In Fig. 14(a) , 

Figure 14(b) shows that the corresponding scan losses of 5.76^ and 3.52% 
respectively. Based on the DOE document (2, p. 27], the antenna/subarray 
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mechanical alignment of HK 3 arc minutes was specified; and in Ref, 3, 
pp, xlv and 1,3 the mechanical pointing and slope accuracy of 2 arc 
minutes and 3 arc minutes (during all phases of operation) were scared, 
respectively. The scan loss corresponding to a slope error of 3 arc 
minuces (or .05") Is 1.79X* Hence the 10 m oscillation Is too severe 
for the antenna structure, 

From Fig. 14(a) p the maximum local slope Is much greater 
than the RMS average; for Instance p It is greater than .13° at t ■■ 10 sec. 
The maximum out-of-plane displacement Is about + 1 m as shown In 
Fig. 14(c). 


It Is Interesting to note that the energy (envelope) of 
the flexible modes leads the energy of the rigid body modes by 1/4 cycle 
of the excitation period as indicated In Fig. 14(d). 

For the purpose of comparison ; two other forcing functions 
of the same amplitude were applied at the antenna Interface. These 

functions are 6832.8 sin (.00675 t) and 6832. 8U(t). 

* 

The responses of cosine and sine excitations are drastically 
different. Figure 15 shows that the latter is extremely smooth compared 
with the former. This is also clear from the energy plot. Fig. 15(d), 
where the flexible modes contribute very little to the system energy. 

In contrast, Fig. 14(d) shows a much higher proportion. Again, the 
dominant modes are 7, 8, and 9. Figure 15(c) shows that the out-of-plane 
displacement has reduced to about half as much and, from Fig. 15(a) and 
(b) , the local slopes and the scan loss are both staying within the 
requirement. 


A step function of the same magnitude, has excited the 
flexible modes in much the same manner as the cosine function has excited 
them. By comparing Fig. 14(b) and Fig. 16(b), for Instance, the two 
responses share approximately the same envelope which is slowly decaying 
with the time constant of the dominant modes, l.e., approximately 1675 
seconds . 
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Table 4, 


The KQsults of these three cases are sumnartxed in 


It is important to point out that in the contest of controlled 
motion! signal shaping is extremely Important and so is the thrust cut-off 
time (refer to Section 3,3,40) , 

B, Surface Response to Collector Bending Disturbance due to Thermal 
Distortion 

In Af the main purposes are to explore the effects of signal 
types on the flatness of the antenna surface. Here! two more cases are 
discussed relative to the orbital thermal effect. It is discussed in 
Section 3.2 that the temperature gradient in the collector can build up 
almost immediately after the SPS moves out from a full shadow and it 
is estimated that for the material with the coefficient of thermal 
expansion (a) of 10” m/m/®C and a temperature gradient of 100®C front-to- 
back! the collector bending can reach the amplitude of 6.67 m. With 
the first collector bending mode dominant, the corresponding force will 
be 

F « 4558.5 cos (.00675 t) N 

_3 

and the corresponding acceleration will have an amplitude of .031 x 10 g. 

The simulated results indicate that this disturbance will 
cause more local surface warping than specified. From Fig. 17 (a) -(d) 
the maximum PIS slope exceeds 0.06“ and the scan loss is about 2.7%. 
Therefore, an ct of 10”^ m/m/“C with 100“C gradient is not acceptable. 
However, an excitation force of 

F - 3417.2 cos (.00675 t) 

corresponding to * 5 m (or 1/4000 of the length of the collector 
structure) will be within the requirement as Indicated in Fig. 18 (a) -(d), 
where the RMS slope error is about .046“ and the maximum scan loss is 
1.45%. The maximum out-of-plan displacement is 0.5 m. From eq. (34) and 
eq. (35), the corresponding a for AT ■ 100“C is .75 x 10”^ m/m/“C, 
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Table h. Comparative Perfomance Analyale. Z •• 10 m 


Applied Force 



6032.8 coa(.00675t) 

6832.8 sln(.00675t) 

6832.8 U(t> 

.0465 X 10"^ 

.0465 X 10"^ 

.0465 X 10”^ 

1000 

1000 

1000 

Modes 7»B,9 

Modes 7,8,9 

Modes 7,8,9 

,091* 

,049“ 

.0905“ 

5,76 

1.73 

5,72 

,069“ 

1 

- 

.0705“ 

3.52% 

- 

3,36% 

1.01 

1 

.56 

1,02 


Table 5. Analysis o£ Performance: Dlstorclona Caused by 

Thermal Disturbance, o ■ 10”^ and .75 x 10“<» m/m/‘'C 


Coefficient of Thermal Expansion 
m/m/®C 

Front-to-Back Temp. Grad., ®C 
Thermal Bending Amplitude, 

Acc. Amplitude, g 
Acc. Force, N 
Local Slope (BMSL max) 

Scan Loss (SLOSS max) , % 

Max. Di8p3,acement (W max) , m 



-6 

10 ” 

.75 X 10 

100“C 

100“C 

6.67 

5.00 " 

,031 X 10“^ 

.023 X 10”^ 

4558.5 cos(.00675t) 

3417.2 cos(.00675t) 

.06 

.046“ 

2.7 

1.45 

,7 

.5 










TheceCore a valuo of a less than ,75 x 10 ^/*C will be acceptable 
ptovldod that the temperature gradient will not exceed 100*0* It seems 
that on a of 10”^ m/m/*C will provide a good safety margin. These 
results ore summarized In Table 5. 

C, Surface Deformation Subject to Collector Torsional Oscillation 

Torsional vibration at the collector boundary can produce 
large torque at the antenna Interface If the amplitude Is 
high. In section 3.2 It is estimated that, for 1® oscillation, the 
torque Introduced to the antenna will be 

T - 6,453 X 10® co8(,0226t) 
with a correspofidlng couple at 1000 m apart of 

F - 6.453 X 10^ co8(.0226t) . 

Figure 19 shows that the resulting surface distortions are within the 
requirement. The maximum RMS slope Is ,05* and the scan loss Is about 
1 . 8 %. 


The dominant modes In this case are modes 10, 12, and 20 as 
Indicated in Fig. 19. 


with 


For the purpose of comparison, another simulation was made 


T " 8.397 X 10® sln(.0226t) 
or the corresponding couple at XOOO m apart, 

F » 8.397 X 10^ sln(.0226t) 

which is equivalent to an oscillation amplitude of 1,3*. A much smuother 
response was obtained as shown in Fig. 20(a) -(d). The results are 
summarized in Table 6. These results again point out the Importance of 
signal shaping. 
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Table 6. Analysis of Performance: Surface Deformation 

Caused by Collector Torsional Vibration 


Collector Vibration 
Amp. at Boundary 

V 

1.3" 

Angular Acc, rad /sec /sec 

,51 x 10"^cos(.0226t) 

.66 X 10"^sln(.0226t) 

Acc. Torque on Antenna, N-m 

6.453 X 10^cos(.0226t) 

8.397 X 10^sln(.0226t) 

Force Couple at Interface, N 

6,453 X 10\os(.0226t) 

8.397 X lo\ln(.0226t) 

Local Slope (RMSL max) 

.05" 

.042" 

Scan Loss (SLOSS max), % 

1.8 

1.2 

Max Displacement 
(W„max) , m 

.7 

.44 


Table 7. Comparative Performance Analysis of Antenna under 
Rectangular-Pulse-Forced Acceleration 


Forcing Function 

F = 4000[U(t)-U(t-500)l 

T « 2xl0^[U(t)-U(t-500)l 

Distributed Force F^,N 

1000 

1000 

Local Slope RMSL max 

.053" 

.017" 

Scan Loss SLOSS max, % 

2.05 

.21 

Max Displacement W max, m 
z 

.6 

.175 


D . Comparative Performance Analyeis of Antenna Under RectaneuIag^Pulse*- 
Forced Acceleration 

Two cases wore simulated, In the first case a rectangular 

pulse force 

F » 4000 tU(t) - U(t-500)1 

was applied at the antenna Interface and In the second case a rectangular 
torque 

T - 2 X 10* tu(n) - U(t-500)1 
was applied co the aiitenna. 

The objectives of this experiment are to make the following 

two points: 

a. The antenna has greater rotational stiffness than 
translational stlffneBs. This observation may be 
verified by comparing the two cases as plotted in 

Fig. 21 and Fig. 22. The force of 4000 N was distributed 
over 6 points (refer to Fig. 13a) and the same force 
was used to produce the torque of 2 x 10* N-m by reversing 
one of the two 3-point sets of distributed force (refer 
to Fig. 13(b)). Note that the force F^ in Fig. 13 is 
.25 F. The translational acceleration has caused a 
RMS slope error of ,053° and .6 m maximum displacement 
whereas the rotational acceleration has caused only 
an .017° maximum RMS slope error and ,175 m maximum 
displacement. Table 7 shows the comparative results. 

b. The cut-off time of the pulse is important. For instance, 
in the case of the translational acceleration, if the 
pulse were cut off at about 26 seconds sooner (or later) , 
the antenna surface would have been flat after the cutoff. 
This is quite cleor from Fig. 21(e) and (d) , since at 
those times the flexible mode energy was at the minimum. 


3*3.5 


Further Discussion 


In Section 3.3.4 the analysis is based on simulation results. 
However, no nuiuttrical Information is by itself complete. In this section, 
we shall complement the numerical results with further analysis. 


Let q be the modal amplitude of a dominant mode, and the 
corresponding differential equation with the excitation of cosine function 
may be written in the following form, 


q 4* 25 wq + w q " h cos m^t 


(39) 


the solution of this equation with zero initial conditions is, 


q(t) 


, 1 2v, ^-ewt 

(m - 0 ) )h e ’ 

o' 




(- costii /l-c t - 




— 5~ — ~ — 7— rr 8ln(o/l-?^t) 


(u)^-w^^)h 


((ja'^-W^‘^)^+(2i:w Wq) 


2Cb) m 

^ (cos w^t + ^ sin u)^t) 

U)^-W 

o 


(40) 


In the context of this application, t; 1 b in the neighborhood of .005 and 
2 2 

<< u) , (40) may be approximated by the following simplified equation 

q(t) » — 5 * {cos m t - cos mt} (41) 

o 

The first term inside the braces represents the steady state solution and 
the second terra contributes only to the transient response. 

2 2 

The maximum value of y(t) is 1.983h/(w -w ) which occurs at 
t “ 26.36 sec. The transient to the steady state envelope ratio is 
about 2 to 1. 

Since the modes are decoupled in our model, (41) may be extended 
to all the flexible modes; hence for aiode k, k » 7 , . . . , 20, 

{cos u) t - e cos u). t) (42) 

o ic 


yk<t) 


h. 


2 2 
m, -m 
k 0 
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From (42), we have the following observations: 


1. The modal response Is proportional to the constant h^^. 

Since hj^ Is the product of the Inverse of the kth element 
of the generalized mass, the mode shape at the point 

where the force Is applied, and the amplitude of the applied 
force, the modal response Is proportional to these 
parameters . 

2 2 

2. For o)j^ » , which Is the case In this application, 

the modal response Is Inversely porportlonal to the square 
of the natural frequency, 

3. Since the surface deformation at a given point Is the sum 

of the products of the local mode shape and the corresponding 
modal response, the above properties apply directly to the 
surface deformation. 

4. Owing to the fact that in this application the local 
slopes are so small that the slope angles are proportional 
to the local deformations : 

a. For a known system response to a specific disturbance 
level, the responses to other levels may be obtained 
by linearly extending the known results. 

b. For known responses to various types of excitations, 
new results may be obtained by taking the linear 
combination of the known results. 

By application of these linear system properties, repeated 
experiments may be avoided and hence a great deal of time 
and effort may be saved. Figures 23, 24, and 25 show 
the performance characteristics of the antenna over wide 
parameter ranges. These figures were obtained by applying 
the above linear properties , 
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Table 8. Key to Variable Names for Figures 14 to 22 


EFLEX 

EL7-EL20 

ERICID 

ESUM 

RMSL 

RMST 

SLOP 

SLRMS 

SLOSS 

T 

WZMAX 

YL7-YL20 


The sum of energy for the first 14 flexible modes, Joules, 

The modal energy for modes 7 to 20, Joules. 

The sum of energy for the 6 rigid body modes. Joules. 

The sum of the flexible and the rigid body modal energy, 
Joules . 

The root-mean-square spatial average of local slope angles, 
degrees . 

The root-mean-square time average of local slope angles, 
degrees . 

The maximum local slope angle at time T, degrees. 

The RMS scan loss (computed from RNSL) , percent. 

The scan loss of the array, percent. 

Simulated time, seconds 

Maximum out-of-plane displacement of the array at time T, 
meters . 

The modal amplitude for modes 7 to 20. 
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Figure 14. Surface Response to Collector Bending Oscillation 
- 10m, F » 6832,8 cos (.00675t) (Cont.). 
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Figure 15. Surface Response to Collector Bending Oscillation: 

Z - 10m, F » 6832.8 sin (.00675t) (Cont.). 
in 
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Figure 17, Surface Response to Collector Thermal 
o » 10-6 m/m/“C, F = 4558.5 cos (.006 












1 Bending Disturbance, 
s (.00675t) (Cont.). 
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KINETIC AND POTENTIAL ENERGY 
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KINETIC AND POTENTIAL ENERGY 



Figure 21. Surface Response to Rectangular Pulse Acceleration - Translation, 
F = 4000 {U(t)-U(t-500)} (Cont.). 
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KINETIC AND POTENTIAL ENERGY 



KINETIC AND POTENTIAL ENERGY 



Figure 22. Surface Response to Rectangular Pulse Acceleration - Rotation 
T * 2 X 10^ {U(t)-U(t-500) } (Cont.). 










AMPLITUDE OF ACCELERATION TORQUE, lO^N-m 
0.0 16.8 33.6 50.4 67.3 84.1 



Figure 24. MPTS Surface Deformation vs. Amplitude of 
Collector Torsional Oscillation 
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MAXIMUM RMS SCAN LOSS 
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APPENDIX A 


Local Slope and Scan Lose Program 


aRUNf/K SLLBBtJaC2LL*S0lARt!>»2U0/9V9V*tlVb/Ab8 .SJW 
dASGfA SPS. 
dOELEUfC 26. 
aCATfP 28. 

8ASG*A 28. 
fiFORtlS SUbl 

SUBROUTINE GRSTR(ViNN) 

parameter NPLOTS8 

REAL Y(20) »YPL0T<NPL0T»2u) 

COHMON/GRSTRI/YPLOT 

C 

C HERE ONLY THE FLtX MODES ARE RECOROLU. TO INCLUDE RIulL oCUY 

C MODESi CHANGE THE RANGE OF K IN THE NEXT INSTRUtMON TO i..;U 

C 

DO 20 K::7t20 
20 YPLOT(NN»K)r Y(K) 

Return 

ENTRY GRR 

call M0UT( YPL0T»ti»8»<;U.3»3H0A = > 

00 30 IdtNPLOT 

URITe(28t40) ( YPLOK IvJliJri >10) 

30 CONTINUE 

DO 35 isifNPEOT 

URITE<2e>40) ( YPL0T(l>J)>Jril>20) 

35 CONTINUE 
AO FORMATi L0E12.6) 

ENOFILE 26 
REWIND 2b 
RETURN 
END 

6F0R>IS SUb2 

SUBROUTINE 6£T0AT(RX>RY#f<2»TXfTYfTZf IPES»lbRm 

Parameter NP=i66>NM=2o>Nt=6i 

REAL RX(NP).RY(NP)>K2(NP) 

REAL TX(NPfNH) f TY(NP>NM) »TZ(NPfNMI 
INTEGER IGRD<NP)fIPt(NE>A)>IPES<Nt>A) 

C 

C NP=NO OF NODES 

C NM=NO OF MOOES 

C NE=NO OF PLATES ELEMENTS 

C 

C RXtRY >R2=C00RD OF GRID POINTS 

C UX>UY>U2=DEF0RMATI0N OF THE GRID POINTS 

C TXtTYf T2 = EIliEN VECTOR MATRICES 

C 16RD=ARRAY OF GRID POINT ID 

C IPE-THE MATRIX OF THE PLmTE ELEMENTS 

C IPES=THE SEO NO ARRAY FOR THt GRID POINTS IN IPE 

C 

C READ GRID POINT FROM SPS.MPTbP 

C 
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Appendix A— Local Slope and Scan Loaa Pcograa 


00 10 1 = 1^4 

RCA0(&»n»£RR:10) RK( 1) «RY (1) tR2(I ) 
to CONTINUE 

00 20 I=1»NP 

RCAOtSfll) IGR0(UfR)(n)»RYU>tRZm 
11 rORnAT(lSt3F15,4l 
20 CONTINUE 

REAP EIGtNVecTORS FROM SPS.MFTEVT 

DO 40 K=lfNM 
00 30 I=lf3 

READ<5.3l.ERR=30) TX ( 1 t K ) »TY U tK) f TZiltKl 

30 CONTINUE 

31 FORMAT(0XtEl2.7f 2E1S*7> 

00 40 IC1*NP 

READ ( 5» 31 1 TX ( 1 *K ) t TY M tK ) »T2 ( I tK ) 

40 CONTINUE 

READ PLATE ELEMENTS »-HOM CAROS 

REa 0( S»51) ( ( IP£(1 *0) »0=1 *4) « icltNE) 

51 FORMAT <161 5) 

test 

uRITE(6fnl) ( ( IP£( I»J) ij=lt4) tl=l«4£ I 

determine SEQ no for IPE 

00 60 IcltNE 
DO 60 KcltNP 
IPESCI »U=IPE( If U 
DO 60 J=2t4 

lF<lPECr*J>.EO.tGRD<K) 1 iPESUtJlcK 
60 CONTINUE 

TEST 


URlTE(ct5U <(1PES(I»J> tJciiA)i 1=1 tKE> 

RETURN 

ENO 

tIS SUb3 

SUBROUTINE MSSLUPES »NE tUX »UY rUZt I CNT »RMST P»C 1 » ABCl t SLOP » 1 SLOP. 
I RMSL.RMST.'ROSS.SLRMS ) 

THIS SUBR COMPUTtS THE LOCAL SLOPES OR ThE DEVIATIONS OF LOCAL 
NORHAT from ITS UNDISTURBED NORMAL VALUE. SINCE THt ANliLES ARl 
small. DOUBLE PRECISION FOR SOME VARIABLES IS REOURED 
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Appendix A-- Local Slope and Scan toaa Program 


IPEScTHE SEQ NO ARRAY FOR T Ht GR 10 POINTS OF THE PLATE ELEMENTS 

NE=N0 OF PLATE ELEMENTS IN ThlS COMPUTATION 

UXfUY»UZ=THE INSTANTANUOUS COORDINATE VECTORS OF THE GRIDS 

SLOPSMAX SLOPE IN THE INTERVAL 

ISLOPsPLATE ELEMENT NO OF MAX SLOPE 

RHSL = SPATIAL AVERACEt RMS SLOPE ERROR. DEGREE 

RHST=RUNNING average. RMS SLOPE ERROR. DtGREE 

DIMENSION IPES( 61 .4 ) . UXd 66 ) .U Y ( 1 66) .UZ ( 1 66 > 

DOUBLE PRECISION A.b.CtCC.ABw 
INTEGER H 

POL=Pl*DL/LAKfaOA*(PI/ieO.) . THETA IN DEGREE 

DE=10.39 METERS ». 

LaMDA = .1224 meter <FRE0 = Z.45 GHZ) 

SLOSS IS THE PERCENT SCAN LOSS OF THE ARRAY 

SLRMS IS the scan LOSS COMPUTED bAStu ON RMS SCAN ANGLE 

PDL24*654375 

SLOSSsO. 

SLOP=0. 

RMSL=0. 

DO 60 1E=1.NE 

H=iPESCIL.2i 

J=IPES(IE.3» 

K = IPES(IE.4) 

MEN! H. J. K ARE THE 1ST. 2ND. 3RD POINT THAT DtFINES PLATE ELEM IE 
A1=UX( J)-UX(HI 
A2=UY (a)-UY(h> 

A3=UZ(J)-UZ(h) 

61=UX(K>-UX (H) 
t)2=UY<K)“UY(h) 

B3=UZ (K >“UZ(H) 

Ar-A3AB2+A2*b3 
6=A3*BI-A1*B3 
C = -A2*Bl-»-AlAb2 
CC=OABS(CI 

ABC=PSQRT< A*A*H*b »C*C> 

COMPUTE LOCAL SLOPE ANGLt IN UEGHLES * 

TEMP=SNGL( (DAC0S(CC/AHC>»*57.2957bDU 
IF ( TEMP. LE. SLOP) GO TO 70 
SL0P=TEMP 
ISLOP=IPESat.l ) 

Cl=SNGL(C) 

ABC1=SNGL( ABC) 

70 CONTINUE 
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Appendix A»- Local Slope end Seen Lose Progren 


C 

c COHPUTt tLEMbNT PATTERN* REClANaULAR APERTURE ANU UNIFORP 

C ILLUNINATION 

C 

ir<TENP.LQ»0) GOTO bO 
AN6=PDL*TEHP 
TENPI=SIN(ANG)/AN6 
SLOSS = SL0SS<^TeMPl*TtPPl 
rhsl=rksutehp*temp 
au continue 

SEOSSs(l.-SLOSS/NE)'‘JOO. 

IFUCNT.EO.O) SLOSSeU. 

RHSL=RKSL/NE 

RMST = RRSTP*RMSTP»ICNT>RMi»L 

COMMENT RMSL IS THE RMS SLOPE - SPATIAL AVERAGE* RMST IS THE TIME AVE 
iCNTPrlCNT^l 
C 

RHST=SQRTtRMST/ICNTP> 

RMSL=SQRT(RMSL) 

C 

C COMPUTE RMS SCAN LOSS 
C 

ANG=P0L*RMSL 

TEMP1-SIN( ANGJ/ANG 

SLRMS=U . 0-TtHP I *T£MP 1 > * 100 . 

RETURN 

END 

OaSGiA csslatran. 

aXQT«G CSSL*TRAN.CSSL 

PROGRAM slope OF lOOO-METER SPS ANTENNA 
INITIAL 

array U<2U>rY< 20> .YU< 20) »YP UO) »YUP <JjO) tZ (20) tUO (20 ) i RH0(20 >*E (20 ) 
array U(20)»bbTA(20) »C8<20|*L1(20) *S1 (20)»TPL0T(ti) 
array RX( 166) *RY( 166) *RZ(i66) tTX (166 t20) *TY ( 166* 20) t T Z (1 66t 20 ) 
array UX( 166) rUY( 166) *UZ(166) *UU(3) *VV(3) *UW(3) 

ARRAY WXU66) *UY(166) *WZ(166) 

ARRAY MS(20)»PF(20) 

integer KK *NN»NPL0T*1P fNEflCNT 

INTEGER IGRDU66) .IPES( 6 1 *4 ) » ISLOP *1 PMAX 

DT=2.5 

TF=750. 

COMHENT THE MAX NO Of PLATt ELLHENTS IS 61 
NE = 61 

COMHENT applied FORCE IN NEWTON 
F0RCE=854.3 
CONSTANT NPLOT=d 

CONSTANT U=0.f0.*0.*0.*0.*0. *.1106toS6 *.1191 V97* .123 192T* .21552SH* ... 
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Appendix A—>LocaI Slope and Sean Loes Program 


.2170462».2!Ve639i>f .27ai*l9t»»»2t«691b5t .29 * .29b0 33d 

• 313267at.3t7&9a4*.32 33IOd*«324a9!>b 
constant nS=3UtO000.»3900U00 ••14600000. *208 0000 *»7lb0000*» ... 

1780000. tlS20000.»lbl0LQ0.t2150UUO,»1340OO0.» ... 
449000U. •3980oat}.t2400000. #7250000. f544a0l)0.» ... 
30000Q0.* 10200000 .tl 120000.* 261 0000. *3260000. 
CONSTANT TPLOT=100.»200.»300.f400«.*500.»60D .*7u0.*600. 

00 L4 KK::lf2b 
Y0<KK )=0. 

V(KK)=0. 

ZfKK)=.00S 

L4. .CONTINUE 

00 Lb KK=1*20 

C8(KK)=SQRT| 1.0 - ZiKKIi^iKKl) 

BETA(KK>= Z(KK)/Cti(KK) 

UO(KK)r U<KK>«C6(KK ) 

C1(KK)= COS(UO(KKIeOTIl 
SKKKlx SIN(U0(KK )a(j7) 

RHOtKK »= EXPI-ZfKK)*y(KK)*uTI 
Lb. .CONTINUE 
NN=1 
SLOPED. 

ISLOPsO 

RNSLrO- 

RHST=6. 

C=0. 

ABC=0. 

ICNTSO 

SL0SS=0. 

SLRRS=0. 

COHHENT CALL GETDAT TO GET DATA PROP SPS FILLS ANU FROP CAPOS 
CALL GETUATIRXfRYfRZtTXtrYtTZt IPESf 1 UROI 
COHNENT initialize OlSPLACtNENT 
DO LIO IPsltlbb 
U)((IP)=0. 

WYOPI^O. 

UZ<IPT=0. 

LIO. .CONTINUE 

DO LiS KK::l*20 

PF(KK) = tTZI136»KK>*TZI53fKK)« .5*(T2a 40»KI(I ♦TZI iSb* KK )♦ 
rZ(58tKM4TZfT7fKKMI/RSIKK ) 

LIS.. CONTINUE 
END 

OYNANIC 

GINTERVAL 01^2.5 

rFIT.GT.TF >60 TO FIN 
IFIDT.NE.CII GO TO FIN 
00 Lib IP=l»lbb 
UX(IP>=0. 

UY(IP»=0. 

UZ(IP>=0. 

Lib*. CONTINUE 

00 Lie KK=lt20 

U(AK>=PFII<l(>*FORCE*COSi .U067&*T> 
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Append.^x A'^Local Slop* and Scan Loa* Progran 


tie. .CONTINUE 

00 L20 KK=lt20 
VP«KK)sY(KK) 

TOP IKK rOIKK ) 
t20.. CONTINUE 

DO t2(< KK=lt20 
iriKK.GT.ei bO TO L23 
CONNENT COHPUTE RlblD BODY NODES 

YIKK)=YPIKK) « YOPiKKMbT ♦ lU IKK » *01 *0T/2 .01 
YOIKK) s YOPIKKI ♦ UIKKUDT 
GO TO L24 

L23.. CONTINUE 

CONNENT CONPUTE FLEXIBLE MOOES 

VCKKI=RHO(KK|a(Cl(KKU BtT A IKK laSU K KH* Yp I KK) ♦ ... 
KH0IKK>*S1<KK>*YLP|KKI/U0IKKI« ... 

U(KK1*U.0-RHC(KK)*IC1 IKKHBETA IKKlaSl IKK ))>/ IW IKK ) •W IKK H 
Y0lKK)=-RH0(KK)«<W(KK)*ISllKKl7CBIKKnaYP IKKI ♦ 

RHOtKK)«ICllKKI>BETA(KKIaSl{KKn* YDPIKKIt ... 
UIKK)*RHOIKK)aSl IKK 1/ IbIKK CB IKK H 

t2«.. CONTINUE 

CONNENT CONPUTE bRIO POINT OISPLACtNENT 
CONNENT 

CONNENT TO INCLUDE RIGIObODY NOTION KK IN US SHOULD START AT 1 
CONNENT •***♦ 

DO L25 KK”?»?0 

DO L25 lP=lfl66 

WXIIP )=UXI IP) + Y(KK >*TXIIP»KK1 

UVilP ) = WYIIP» + YIKK)*TYIIP*KR> 

UZ(lP)=UZ(IP)^y|KK»*TZIIPfKK> 

L2S. .CONTINUE 

CONNENT GRID POINT POSITION 

UNAX=0. 

TENPSO. 

IPNKX=0 

DO L26 tP=l»l66 
UX(IP)=RXI1P)-<'UXMP) 

UYIIP>=RYIIP)-»UVIIP) 

UZIiP)=RZaP)*UZ( IP> 

CONNENT CONPUTE NAX OUT-OF-PLAN OSSPLACEHENT 
» ABSUZcABSIUZIIP) ) 

IFCTENP.GE. AbSUZI 60 TO L26 
MZNAX=UZI1P) 

TENP=A6SIU2HAX) 
rPNAX=IGRDIlPi 
L26. .CONTINUE 

CONNENT CONPUTE SPATIAL AVE RNS SLOPE ERROR ANu RUNNING AVE RMS SLOPE cFROR 
RNSTPsRMST 

call MSStI IPESiNEtUX»JYtUZ»lLNT vRNSTP »CfABCtSLOP»ISLOPtRNSLt ... 

RHST.SLOSS.SLRMSI 

ICNTsICNT^l 

IFINN.GT.NRLOTI 60 TO L3U 
IFIT.LT.TPLOTINNH GO TO L30 
CALL GRSTRlYtNNI 
NN=NN+1 
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Appendix A--I«oc«l Slope and Scan Loaa Prograx 


L30..C0NTINU£ 

CRIGIOSO, 

CFLEXsO. 

00 L50 KKsl^;20 

EfKKI s ¥0(KK)*VO(KK) « UfKK ) aU< KKI * VI KK I eV(KK ) 
EIKK>=.5*HS(KK)ae(KK) 

L5U.. CONTINUE 

00 L55 KKslf6 
£RIGI0=ERIGIU«EIKK» 

L5S.. CONTINUE 

00 160 KKsTfPO 
EFLEXrEFUX^lIKKI 
160.. CONTINUE 

ESUN=ERIGIOHFL£X 

VU=VI|> S YLisY<2) S yU = Y|3» % yimUHi % YL5 = V(5I 

Vl6=YI6> I Y17=Y(7) i YU>=YU» f VL9=Y<9) i YUOsVIlO) 

VU1=YIIU • YU2=YU2) * YU3=YU3) I YU4=Y(U) * YUSsYUSl 

VU6=Y(16) » YU7=YU7) I YLl8rV|ie» S YU9=Yll9> t YL20=YI20| 

ELltEll) $ Ei.2sE(2) S EU=E(3) ft U«~£(4) » EL5=C(S> 

EL6SEI6) ft EL7=Ern ft tUd=E(tt» ft EL9=E(9) ft ELlOsEllOl 

ELllsEIll) ft EU2-EU2J ft EU3 = £U3) ft EU4=EI19) ft fcL15=£(l5) 

EU6=eil6> ft EU7=£(m ft ELiarEllO) ft EU9=ta9> ft tL20=tl20) 

OUTPUT YLlfYL2»YL3»YL4»YLSiyL6fYL 7i YL 6 fYLY »yUOfYUl #yLl2*yU3* ... 

VU<l»VLt5»VLl6»YU7»YUe»YU9» YL.!OtfiFSTtR«SUilSLQP»SLOP» . .. 
C»ABC»SL0 5SfSLRHS*lPHAXfW2HAXtESURfEKlGIUiEFLEX 
PREPAR YU*YL2*VL3»YL4»YL5»Yi6»YL H YL8 .YL9 »YLlO*YLll #YU2* YU3» ... 
YU4»VUS.YU6fYLl7tYLldfYLl9.YL20tRXSTtRHSLtSL0p*LSURf ... 
ELltEL2tEL3»LL4.a!>.bL6*EL7fU6»LL9*EU0*El ll tEL12»EL13t ... 
£U4tEL15»EU6«EU7fEL18*CL19*a2i)»Sl0SStSLRHSfUZNAXt ... 
ESUH.ERl610»tfLEX 
OERIVATIVE GRP 
YOt=i.O 

YIsINTeGIYuI tC.G) 

ENO 

END 

terhinal 

FIN..CONTINUE 
CAU 6RR 

ENO 

ENO 


■Vi 
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rACT0R«*75> 

LABCW nOOAL COOROIRATCS 

GRAPH T.YUTfYtA.YLV 

UBCL HOOAL COOflOlRATCS 

GRAPH T*YL10tYUl>YU2 

label modal COOROIMATCS 

GRAPH T>YLL}«YLIR*YLlS»yLlft 

LABEL MODAL COORDINATES 

GRAPH TtYH7*YLIB*YL19.YL23 

label MNETIC AND POTENTIAL ENERGY 

GRAPH T#EL7«ELB»CL9*EL10 

LABEL KINETIC AND POfENTlAL ENERGY 

GRAPH TtELI liCL\2*£Lt3«CLU 

LABEL KINETIC AND POTENTIAL ENERGY 

GRAPH TfSLl5f£Ll6fELt7»CLtd 

LABEL KINETIC AND POTENTIAL ENERGY 

GRAPH T*eLI9»EL20 

label kinetic and POTENTIAL ENERGY 

graph TtESUMfENlGlOtEFLEX 

LABEL MAX OUT-OP -PLAN DISPLACEMENT 

graph r.RZMAX 

LABEL ANG DEV IN DEG - RMS RUNNING AVEt R^G SPaTIAL AVtt •’AY 
GRAPH TtRHSTfRMSL.SLOP 

label PERCENT SCAN LOSS AND RMS SCAN .OSS 

GRAPH T»SLOSS«SLRHS 

STOP 

BEOiV 2B.tSPS.2B 
BE OP 

BPIN Mri 
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Appendix B 

3-D Plotting Program of Spatial Surfaces 


aRUN«/R SPSAN»J6C2GLtS0LARf3f200/9999f «t93/65B .SJU 

aASGfA SPS. 
arORtlS MAIN 

PARAMETeR NP=166tNH=20tNC=272tNPL0T- B« NM2=1Q*NPL0T2:1& 
real RX(NP >tRYfNP)»RZfNP) tV(NM>fyE0U(NPLUT2*NM2) 

REAL UX(NP) *UY(NP)«UZ(NPtfTt<NPtNMMT2(NPtNMI»T3<NP«NM) 
INTEGER S£0(NCltS£01(NPI 
C 

C NP=NUHB£R OF NODES 

C NM=NUMBER OF MODES 

C NCcNUMGER Of plot commands 

C NPL0T=NUMBER OF PLOTS 

C NM2=NM/2 

C NPL0T2=NPL0T<*2 

NMZizNMZ*! 

0R=3,t416/iaQ.Q 

call plots 

CALL PL0T<6.a»A.0*-3» 

CALL FACTOR (0,51 

c read grid 

DO 5 I=lf4 

READ(5t333feRR=5) RXf I) f RY( I) « RZ ( II 
5 CONTINUE 

DO to 1=1 fNP 

REA0(5f33J) SEQI(n»RX<nfRYUI«RZm 
333 FORMATtl5»3F15.U 

10 CONTINUE 

C READ PLOT COMMAND SEQUENCE 

READ(5>;334>(SE0(J)fJ = l tNC) 

339 FORMAT (1615 I 

C READ EIGENVECTORS 

DO 20 K=1»NM 
DO 22 1=1*3 

READ (5*335* ERR = 22i T1 ( I *K ) *T2 ( I* K) * T 3( I* K I 
22 CONTINUE 

335 FOR M AT ( ax * E 1 2 . 1 * 2£1 1 ) 

DO 20 I=1*NP 

20 READ(5»335> T1 ( I »K I *T2 ( I*K I • T3 ( ItK 1 

REaO( 5*336) ((Y£QU( l*J}*J=t*NM2l *1=1 *NPLoT2) 


336 F0RMAT(10C12.6I 
DO 30 I=1*NPL0T 


DO 25 J=1*NM2 
Y(J)=YEQU(Itsl) 
25 CONTINUE 


DO 27 J=NM21*NM 
Y( J)=Y£QU(NPL0T»I*d-NM2l 
27 continue 


C 


Appendix B— 3-D Plotting Program of Spatial Surfaces 


CALL 6RSUN( YfRXtRYtRZfn iTZtTltUXtJY tUZ) 

Call oRAUfuxtUVtUZtSEatSCQii 
CALL PLOTaa.OfO,Of-3l 
30 CONTtNUe 
C 

CALL PLOTUO.OvO.Of 999) 

50 CONTINUE 
STOP 
END 

aFORfIS SUHl 

SUBROUTINE GRSUN) ViRX»RY»RZ»Tl.Ti»T3»UXtUYfU2) 

PARANETER NH=20»NP=I66 

REAL TlINPfNN) tTZINPtNN) fT3(NP«NH) 

real yiNM>tRX<NP)rRV(NplfRZINP) »UX (NP ItUVlNP) •UZ<NP) 

SCL1=0.005 

SCL2=2. 

DO 5 IrltNP 
UXII)=0.0 
UV<I)=0.0 
5 uzm=o.Q 
00 10 K=:ifNR 
00 10 I=ltNP 

uxm= UXU)-» Y(K)*T1( I.K)*Sa 2 
UYa)= UYtll* Y(K)*T2( I.K)*SCL2 
10 uzci)= UZ(I)^ Y(K)*T3{ lfK)*SCL 2 
LO 20 IrlfNP 

UX(I)= RXII)«>SCL1 * UX(I) 

UY<n= RY(I)*SCLl * UYin 
20 uzin= Rzm*scLi ♦ uz(i) 
return 

END 

aroRvis SUB 2 

SUBROUTINE DRAU( UXiUYfUZ tSEQrSEQ U 
PARAHETER NP= 1 66 . NC=272 
REAL UX(NP),UY(NP).U2(NP) 
integer SEQ(NC) f FLAG*5EQKNP) 

00 10 1=1 »NC 

J=1ABS(S£G(I)> 

00 5 K=1*NP 
IFfSEQKKI.EQ. J) KK=K 
5 CONTINUE 
FLAG=2 

IFISEOa ).LT.ai FLAG = 3 
X = UX(KK) 

V=UY1KK) 

Z=UZ(KK) 

call TRANS<Xf YtZ»XP»VP) 

CALL PLOT (XPfYPrFLAGI 
10 CONTINUE 
RETURN 
END 


Appandlx U--3-D Plotting Ptogtam of Spiitlal Surfaced 
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